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ABSTRACT 

We present our analysis of the X-ray variability of two ultraluminous X-ray sources (ULXs) 
based on multiple XMM-Newton observations. We show the linear rms-flux relation is present 
in eight observations of NGC 5408 X-1 and also in three observations of NGC 6946 X-1, but 
data from other ULXs are generally not sufficient to constrain any rms-flux relation. The 
presence of this relation was previously reported in only two observations of NGC 5408 X-1; 
our results show this is a persistent property of the variability of NGC 5408 X-1 and extends to 
at least one other variable ULX. We speculate this is a ubiquitous property of ULX variability, 
as it is for X-ray variability in other luminous accreting sources. We also recover the time 
delay between hard and soft bands in NGC 5408 X-1, with the soft band (<1 keV) delayed 
with respect to the hard band (>1 keV) by up to ~10 s (~0.2 rad) at frequencies above ~few 
mHz. For the first time, we extend the lag analysis to lower frequencies and find some evidence 
for a reversal of the lag, a hard lag of ~1 ks at frequencies of ~0.1 mHz. Our energy-resolved 
analysis shows the time delays are energy dependent. We argue that the lag is unlikely to be 
a result of reflection from an accretion disc (‘reverberation’) based on the lack of reflection 
features in the spectra, and the large size of the reflector inferred from the magnitude of the 
lag. We also argue that associating the soft lag with a quasi-periodic oscillation (QPO) in 
these ULXs - and drawing an analogy between soft lags in ULXs and soft lags seen in some 
low-frequency QPOs of Galactic X-ray binaries - is premature. 

Key words: X-rays: general - X-rays: individual: NGC 5408 X-1 - X-rays: individual: 
NGC 6946 X-1 


1 INTRODUCTION 

In recent years substantial progress has been made in understanding 
the nature of ultraluminous X-ray sources (ULXs). These objects 
are X-ray sources located within, but displaced from the nucleus 
of, nearby galaxies, that d isplay observed X-ra y luminosities in ex¬ 
cess of 10^^ erg s“’ (see lpeng & SorialUollL for the most recent 
review). In particular, there is strong evidence for three relatively 
nearby ULXs containing stellar-mass blac k holes (Mbh ~ 1041a) , 
and so accreting at super-Eddingt on rates jMiddleton et aljr2013l : 
IlIu et alj2013l : llVtotch et alj|2014h . The last of these objects is par¬ 
ticularly important as it directly links super-Eddington emission to 
the peculiar X-ray spectrum displayed by many bright ULXs (see 
e.g.. IStobbart et aT]|2Q0d : iBachetti et alj|2013h . and so supports 
the notion of som e ULXs accreting in a new, super-Eddingt on ul¬ 
traluminous state ([Gladstone et alj|2009l : ISutton et al.ll2013l) . The 
physics of this state are only just emerging, but it appears that the 


* E-mail: lorena@iaa.es 


variety of X-ray spectra and the coarse variability properties of in¬ 
dividual objects may be consistent with a model in which massive, 
radiatively-driven winds are launched from a geometrically thick 
accretion disc, bloated by the advection of its hot inner regions to¬ 
wards the black hole, and a combination of the collimation of the X- 
radiation emerging from the innermost regions of the accretion flow 
by the optically thic k wind and the viewing angle of the observer 
dictates what is see n jPout anen et al.l2007l : lKinj20()^ : ISutton et al.l 


I 2 OI 3 I : [Middleton et ^ 1201^ However, this may not explain all 

ULXs, which we now are certain to be a heterogeneous population 
after the discovery of pulsations from a lu minous ULX in M82 , 
demonstrating that it hosts a neutron star dBachetti et alj|20l3) . 
Additionally, there is evidence that at least some of the brightest 
ULXs may still harbour the long- sought intermediate-mass black 
holes (Mbh ~ 100 - IOOOOMq, e.g. lFawell et al.l2009l : ISutton et al.l 
l2012l : l^ezcua et al.l2015h . Much work therefore remains to be done 
to understand both the composition and the accretion physics of the 
ULX population. 

X-ray variability can be a powerful tool for investigating lu- 
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Figure 1. rms spectrum (black circles) and mean spectrum (blue squares) of (left): NGC 5408 X-1 using the 2010 and 2011 observations, and (right): NGC 6946 
X-1 using the 2007 and 2012 observations. The lower panels show the ratio between the rms and the mean spectrum, i.e., the fractional rms spectrum. We used 
a time bin size of 10.4 s respecting the “good” time interval list, and segments of 10 ksec duration for the estimation of the rms (see text). 


minous^_accasting_blac^ioile_bitiaries_(BH ^Reinillard & McClin - 


tock 20061) or active galactic nuclei (AGN. Vaughan et al.ll2003bh . 

Studies of X-ray variability have been undertaken for a few ULXs 
using differentag£roacheSj_^om_changesjndheirlightcurves_(Sut- 
ton et al.l2013[). the study of q u asi-periodic oscillations (QPO , e.g., 
IPheerai & Strohmaveilfior^ ICaballero-Garcia et al.ll2013ll . the 
rms-flux rePtion or time delays between different energy bands 
jHeil et al .11 20091 IPe Marco et alj|2013bh . However, the study of 
variability in ULXs is hampered by their relatively low count rates, 
and usually requires long exposures in order to recover the statis¬ 
tical properties of their variability. For example, the rms-flux re¬ 
lation or time la gs have to-date be e n dete cted in only one ULX, 
NGC5408 X-l. lHeil & VaughanI l l2010t) first showed that this 
source follows the rms-flux relation from two XMM-Newton ob¬ 
servations. This property of the X-ray variability seems to be ubiq¬ 
uitous among Galactic black hole X-ray binaries (XRB) and als o 
AGN dUttlev & McHardvl200ll:rUttlev et alJ200^:lHeil et al.l2012h . 


The discovery in NGC 5408 X-1 demonstrates a strong connection 
between at least this ULX and th e better-understood accr eting black 
hole systems in XRBs and AGN. lHeil & Vaught boioh also stud¬ 
ied time delays between the soft and hard X-ray energy bands in 
NGC 5408 X-1 and found a soft lag, i.e., variations in the soft pho¬ 
tons lag those in h ard photons at mHz frequ encies. This result was 
later confirmed bv IPe Marco et alj j2013bh . who used six XMM- 
Newton observations to study the QPO and the soft time lag. 


Here we revisit the XMM-Newton observations of ULXs, con¬ 
centrating on NGC 5408 X-1 and NGC 6946 X-l. Our focus is on 
the frequency dependent time lags (extending to lower frequencies 
than previously published) and confirming and extending the one 
published example of the rms-flux relation in a ULX. This paper is 
organised as follows: in Sect. [2] we describe the data used for the 
analysis and the data reduction, in Sect. we present the rms-flux 
relation for the two ULXs, and in Sect.|4]the cross-spectral analy¬ 
sis and the results of the coherence, time and phase delay, and the 
phase lag spectrum. Finally, the results derived from this study are 
discussed in Sect.O 


Table 1. Observational details. 


Object ObsID Date T“ T* 

(ksec) (ksec) 


NGC 5408 X-1 


NGC 6946 X-1 


0302900101 

0500750101 

0653380201 

0653380301 

0653380401 

0653380501 

0723130301 

0723130401 

0500730201 

0500730101 

0691570101 


2006- 01-13 
2008-01-13 
2010-07-17 

2010- 07-19 

2011- 01-26 
2011-01-28 
2014-02-11 
2014-02-13 

2007- 11-02 
2007-11-08 
2012 - 10-21 


130 

99 

113 

47 

104 

60 

128 

111 

119 

90 

124 

95 

35 

34 

33 

32 

28 

25 

33 

31 

114 

99 


(a) Exposure time before the background filtering, and (b) duration after the 
background filtering (see text). 


2 DATA REDUCTION 

In this paper we utilise the multiple XMM-Newton observations of 
NGC5408 X-1 and NGC6946 X-L These are among the bright¬ 
est and most variable known ULXs and have some of the longest 
publically available XMM-Newton observations. We have in fact 
performe d much of our analy sis on all 20 ULXs in the sample dis¬ 
cussed bv ISutton et al.l(l20L3h . but for the other 18 sources the data 
were not sufficient to obtain meaningful constraints from the time 
lag and rms-flux analyses. We used XMM-Newton archival data up 
to March 2015. The log of the observations used in our analyses is 
given in Table [T] 

For these sources we rejected the observations with exposure 
times shorter than 10 ksec because the analysis was performed us¬ 
ing continuous segments of this length in order to reach low fre¬ 
quencies. To obtain a high signal-to-noise ratio (S/N) we combined 
data from EPIC pn, MOSl, and MOS2 detectors. We forced MOS 
and pn light curves to have the same tstart and tstop. The event 
files for the source were extracted from circular regions (aperture 
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NGC 5408 X-1 (2010a,b + 2011a,b) NGC 6946 X-1 (2012) 




Figure 2. rms-flux relation in the soft (0.5-0.9 keV) and hard (1.2-7 keV) energy bands for (left): the four observations from 2010 and 2011 of NGC 5408 X-1, 
and (right): the 2012 observation of NGC 6946 X-1. The rms is measured over the 4-50 mHz frequency range from segments of length 150 s (NGC 5408 X-1) 
and 250 s (NGC 6946 X-1). The dashed lines show the 95% “conhdence bands” around the best-fitting linear model. 


radius of 40" and 25" for NGC 5408 X-1 and NGC 6946 X- iQ) and 
for the background from rectangular regions, free of other sources 
and located in the same chip of the source, using the Science Anal¬ 
ysis Software (SA®, version 13.0.0. Standard event patterns (pat¬ 
tern ^ 4 for the pn detector, pattern ^12 for the MOS detectors) 
and filter (flag = 0) were used. 

The subsequent analysis was carried out using IDL softwar^. 
The ULXs are relatively weak sources (in some cases <0.1 ct/s 
in the 2-10 keV band), and so may be overwhelmed by the back¬ 
ground during periods of high and flaring background. Strong back¬ 
ground flares may introduce spurious variability and time lags into 
the data (the spectrum of a background flare evolves with time) if 
not properly excluded from the data. In order to mitigate against 
this, we carefully filtered each observation for background flares 
as follows. A light curve of the background was extracted in the 
1-12 keV energy band using 2.6 s bins and smoothed using a 500 s 
width boxcar filter to improve the S/N. We manually selected a suit¬ 
able background threshold for each observation, above which data 
were considered to be affected by background flares, and excluded 
from further analysis. The quiescent (non-flaring) background lev¬ 
els were slightly different between each observation. We therefore 
chose a slightly different background threshold levels for each ob¬ 
servation in order to leave periods of quiescent background but 
effectively remove flares. The background thresholds were in the 
range 0.028-0.036 ct/s for the NGC 5408 X-1 observations, and in 
the range 0.010-0.015 ct/s for the NGC 6946 X-1 observations. 

Col. 5 in Table [T] shows the total amount of “good” exposure 
time after these high background periods have been excluded. We 
note that the rms-flux analysis (Sect.[3]l and the cross spectral anal¬ 
ysis (Sect. 0 were each carried out using equal length segments 
of continuous good data, but the lengths of the segments used for 
each analysis were different. This is because the rms-flux analysis 
concentrates on the higher frequency (i.e., shorter timescale) varia¬ 
tions (the “red noise” part of the PSD) while for the cross spectrum 
analysis we are interested in the variability properties to lower fre- 


' PN data of obsID 0200670301 and 0200670401 were extracted from el¬ 
liptical regions to avoid the detector chip gaps. 

^ http://xmm.esac.esa.int/sas/ 

^ Available from http://www.star.le.ac.uk/sav2/idl.html 


quencies (i.e., longer timescales). Requiring fixed length segments 
of uninterrupted good time means some small intervals of good 
time were not used; the amount of “good” data used for each anal¬ 
ysis is therefore slightly different (lower) for each observation than 
the “good” duration given in Table[T] 

Finally, we obtained an rms spectrum (defined as the square 
root of the normalised excess variance; see I Vaughan et alJl2Q03bL 
for details) of the source (see Fig.[T](, and checked at which energies 
variability is found in ord er to select the energy bands. The spec¬ 
tral analysis performed bv iMiddleton et al] ^201 1[) showed the need 
for two spectral components to fit the spectrum of NGC 5408 X-1, 
which are separated at ~ 1 keV. Therefore, we will take 1 keV as 
the separation point for the energy band selection. From this analy¬ 
sis we selected the soft and hard energy bands as 0.5-0.9 keV and 
1.2-7.0 keV as these were the bands where the strongest variations 
were found. 


3 RMS-FLUX RELATION 

In this section we describe our analysis of the rms-flux relation in 
these two UL2^SjThe^n^^sisJollGws_clGsel^_die_anal^sis_o^Heil 
& Vaughan dTOlOh and we ref er the interested reader to this paper 
(and that of IHeil et alj|201lh for specific details. Briefly, we di¬ 
vided each observation into continuous 150 s and 250 s segments 
of “good” data, for NGC 5408 X-1 and NGC 6946 X-1 respectively, 
and from each segment computed a periodogram and the mean 
count rate. We needed slightly longer segments for NGC 6946 X-1 
in order to achieve a good rms estimate, since the rms is lower in 
this source (see Fig. [TJ. The periodograms were then averaged in 
eight groups according to the mean count rate, and the rms esti¬ 
mated from each. The rms is the square root of the Poisson noise- 
subtracted variance, itself computed by integrating the average pe¬ 
riodograms over the ■'-4-50 mHz range, using a 10.4 s (four times 
the MOS frame time) time resolution of the light curves. The anal¬ 
ysis was performed for a soft and a hard band. Linear models of the 
form cr = k((F)-C), where F is the flux and k (i.e., dcrld{F}) and C 
are constants, were fitted to the rms-flux data using weighted least 
squares (min x^)- We also estimated the 95% “confidence bands” 
around the best-fitting linear model. We randomly generated 500 
models, each one from the distribution of parameters (specified by 
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the best-fit values and covariance matrix), and extracted the 2.5% 
and 97.5% y values at each x value, and within these we estimated 
the 95% confidence band. These illustrate the uncertainty on the 
fitted_models^^^^_ _ 

Heil & Vaughan ( |201(]|) demonstrated a positive, linear rms- 
flux relation for the harder band data during the 2006 and 2008 ob¬ 
servations of NGC 5408 X-1. Appendix lAl shows the rms-flux rela¬ 
tions for the 2010, 2011 and 2014 observations (the best fitting lin¬ 
ear model parameters were consistent between the two 2010 obser¬ 
vations, between the two 2011 observations, and between the two 
2014 observations, and so the closely spaced pairs of observations 
were combi ned). In agreement with the results for 2006 and 2008 
obtained bv iHeil & Vaughaij fcOlOh . the 2010, 2011 and 2014 ob¬ 
servations also show the rms-flux relation in the harder band data, 
that are well fitted by a linear model (p > 0.05 in all cases), although 
the parameters are slightly different between the three years. These 
detections (in the sense that the gradient k > 0) are significant at 
the 4.8(T, 7.8(T and 6.6(T levels, respectively. We notice that the 
gradients and intercepts between the 2010 and 2011 observations 
are consistent within the 2a level, so we have combined the four 
observations. Fig.|2](left) shows the rms-flux relation for these ob¬ 
servations of NGC 5408 X-1, whose detection is significant at the 
1.2a level in the hard energy band. 

The soft band appears significantly different (see left panel in 
Fig .|2] and Appendix|3; the gradient k is much smaller. The gradi¬ 
ent k gives the fractional rms after subtracting any constant flux or 
rms components; the much lower gradient in the soft band indicates 
that, even after removing any constant component, the fractional 
amplitude of the variability is greatly suppressed. A linear model 
fits the data well (p > 0.05), although all the detections are below 
the 3a level, so the rms-flux relation cannot be confirmed in the 
soft energy band of NGC 5408 X-1. 

We also show the rms-flux relation for NGC 6946 X-1 for the 
first time (Fig.|2] right panel for the 2012 observation). Data from 
2007 were used together because the linear fit was consistent be¬ 
tween the observations (see Appendix lAt. and the 2012 data was 
treated separately. The rms-flux relation is clearly detected in the 
hard band at the 4.5cr and 4.4cr level for the 2007 and 2012 data, 
respectively. In the soft band the detections are below 3a, although 
again we obtain p > 0.05 in every fit with a linear model. 

Therefore, the rms-flux relation is clearly detected in the hard 
energy band of both ULXs, but we cannot claim this relation for 
the soft energy band in these ULXs. 



0.1 1.0 10.0 
Frequency (mHz) 


Figure 3. Results of the cross-spectral analysis between the soft (0.5-0.9 
keV) and hard (1.2-7.0 keV) energy bands for NGC5408 X-1 showing, 
from upper to lower panels: the coherence; phase-lag; and time lag between 
the energy bands (the inset panel shows a zoom-in of the soft lag). The 
phase lag above ~1 mHz is is significantly negative, at around -0.2 rad, cor¬ 
responding to a time lag of ~10 s at ~4 mHz. At lower frequencies the lag 
becomes consistent with zero or even positive lag in the lowest frequency 
band. The middle panel shows a constant model (dot-dashed line) and a 
power law plus a constant model (dashed line) fitted to the phase lag (see 
text). 


4 CROSS-SPECTRUM 


Flere we report our cross-spectral analys is, which was carrie d out 
following standard procedures (see e.g., iNowak et alJI 19991) . We 
extracted the background-subtracted source light curves in the soft 
(0.5-0.9 keV) and hard (1.2-7.0 keV) energy bands using time bins 
of 10.4 s respecting the “good” time interval list. We used segments 
of 10 ks duration, and from these we estimated the power spec¬ 
tra (PSD), noise-subtracted coherence, phase and time lags in the 
Fourierfte 2 uenc^_domaininfiiestandard_wa^_(seeej^Vaughan & 
Nowak' Il997l ; Nowak et al.lll999l ; IVaughan et alJl200^ . The frac¬ 
tional normalisation was adopted and we rebinned over logarithmic 
frequency bins, each spanning a factor ~ 1.7 in frequency. 

A soft or negative lag (i.e., the soft band photons lag those 
in the harf^and}_m_^GC5408_XH^^as_^rst_^egorte^^^FIeil & 
Vaughan ( l2010li from the XMM-New ton data from 2006, and later 
confirmed by De Marco et alj ( l2013bl) using six XMM-Newton ob¬ 


servations between 2006 and 2011. IPe Marco et al.l ll2013bh com¬ 
bined data from the same year as stationarity was observed. Since 
the PSDs_differbetweenJhe^earSj^s^revmusl^^oted_^De Marco 
et al. (l2013bli . we first undertook the analysis separately for each 
year, as in Sect. [3 From integrating the PSD, we found that the 
variability is stronger in 2010 and 2011 (four observations in to¬ 
tal) and the phase lag spectra appear consistent, and thus we com¬ 
bined these data in order to obtain a higher S/N. The results of the 
cross spectrum for NGC 5408 X-1 are shown in Fig. including 
the noise-subtracted coherence, phase, and time lag between the en¬ 
ergy bands. The soft lag is detected at the A.2a level. In Fig.[3]it can 
be seen that at lower frequencies the sign of the time lag changes 
and becomes positive. However, we cannot claim a hard or positive 
lag (i.e., the soft band photons lead those in the hard band) because 
this measurement is at the 2.9cr level. The strength of the estimated 
lag also varies with the choice of energy bands. 

Although a similar time lag behaviour is observed in 
NGC 6946 X-1 (see Appendix|0, the S/N of the data is not enough 
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Observed Energy (keV) 


Figure 4. Phase lag spectrum of NGC 5408 X-1 at low (black; 0.1-0.3 mHz) 
and high (red; 1-20 mHz) frequencies. These show the average lag between 
a given band and the 1-1.5 keV reference band. The high frequency lag- 
energy spectrum decreases with energy, indicating the phase lag becomes 
more negative for widely separated energy bands. The y-axis are different 
for low (left) and high (right) frequencies. 


to confirm the lags, whose detections are significant at the 2.3cr 
(soft lag) and 1.3 it (hard lag) level at high and low frequencies, 
respectively. 

We fitted two different models to the phase lag (see middle 
panels in Fig.j^and Appendix 1^. The first one is a constant model 
(dot-dashed line), and the second one consists of a power law plus 
a constant (dashed line) as phase lag A0 = A/““ -l- B, where A and 
B are constants. The second model fits both sources better, with 
A = 0.05±0.07, B = -0.21 ±0.04, a = \3 ±0.1 Id.o.f = 9.2/7 
and p = 0.24 for NGC 5408 X-1, andA = 0.7 ±0.7, B = -0.6 ±0.6, 
Q- = 0.4 ± 0.4, x^ld.o.f = 1.9/7, and p = 0.97 for NGC 6946 
X-1. With this model the hard lag at low frequencies has a A(/> oc 
/“' dependence, and a constant phase soft lag of -0.2 rad at high 
frequencies. Therefore the phase lag is not constant with frequency, 
nor is the time lag. 

For NGC 5408 X-1 we also computed the phase lag vs. en¬ 
ergy spectrum from the cross-spectra, averaging over the frequency 
bands. Fig.|4]shows the phase lag spectrum from the 2010 and 2011 
observations of NGC 5408 X-1 at low (0.1-0.3 mHz, in black) and 
high (1-20 mHz, in red) frequencies, corresponding to the hard and 
soft lags for NGC 5408 X-1. The 1.0-1.5 keV energy band was used 
as the reference band because this is where the signal is highest, 
i.e, the rms is high. The lags have not been shifted, so zero-lag 
means that there is no time delay between that bin and the refer¬ 
ence band. Similarly, a negative lag means that the bin leads the 
reference band, and positive lags mean that the bin lags behind the 
reference band. At low frequencies (i.e., the hard lag) the phase 
lag increases with energy, from —0.6 at 0.8 keV to ~ 0.3 at 3 
keV. Conversely, at high frequencies (i.e., the soft lag) the phase 
lag spectrum decreases with energy, from ~ 0.2 at 0.8 keV down 
to —0.1 at 3 keV. Note that the y-axis is different for low (left y- 
axis) and hi gh (right y-axis) frequenc ies. The same behaviour was 
reported by IPe Marco et alJ (12013bl) for the time lag vs. energy 
spectrum of the soft lag. 


5 DISCUSSION 

We have shown that a linear rms-flux relation is present in the hard 
energy band of the two ULXs studied here, NGC 5408 X-1 and 
NGC 6946 X-1, and that it is persistent across multiple observa¬ 
tions spanning ~8 years. The limitations of the soft band data are 
such that we are unable to place any interesting constraints on the 
rms-flux relation in the soft band (this is due to the relatively low 
total count rate of these sources, < 1 ct s“*, and the low fractional 
rms amplitude below 1 keV, which means the intrinsic rms is low, 
the Poisson noise is relatively high and the flux range is small). We 
also examined the frequency-dependent X-ray time lags, extending 
the analysis to lower frequencies than previous studies. We find the 
coherence between soft (<0.9 keV) and hard (>1.2 keV) bands is 
high, consistent with 1 at all frequencies. The soft phase lag can 
be recovered for most observations of NGC 5408 X-1 at frequen¬ 
cies above ~few mHz, but the phase lag is not constant down to the 
lowest frequencies we probe (~0.1 mHz) and becomes less nega¬ 
tive at lower frequencies. We find tentative evidence for a hard lag 
(positive phase lag) at the lowest frequencies. 


5.1 The rms-flux relation 


Heil & Vaughan ( 1201 (j) used the 2006 and 2008 observations to 
show that NGC 5408 X-1 exhibits a linear rms-flux relation. We 
show the relation persists through further observations in 2010, 
2011 and 2014. The parameters of the best-fitting linear models 
differ between observations, presumably due to modest changes in 
mean flux and PSD shape. We also detected a linear rms-flux rela¬ 
tion in the 2012 observation of NGC 6946 X-1, demonstrating that 
NGC 5408 X-1 is not unique among ULXs. We did ex amine the X- 
ray va riability of the other ULXs from the sample of ISutton et al.l 
(I 2 OI 3 I) but we were unable to obtain useful constraints on the possi¬ 
ble rms-flux relations. Given that we did detect the linear rms-flux 
relation in the two best observed, highly variable ULXs, we specu¬ 
late that the rms-flux relation is common to variable ULXs but that 
in most other cases the low count rates and shorter observations 
preclude a good detection. 

The rms-flux relation appears to be ubiquitous in luminous, 
accreting_obiectSj_juch_as_jctivg_galacticjuclei_X^GMi, Vaughan 
et al. l2003allbl) . Galact ic X-ray binaries (XRB) w ith black hole or 
neutron star accretors dUttlev & McHydv 200 it Uttle^ 20041) and 
catac lysmic variables dScaringi et al. 20121 ; ~n de Sande et al.l 
l2015t) . The simplest conclusion is that a common physical mech¬ 
anism explains the strong, aperiodic variability, following a lin¬ 
ear rms-flux relation, in all these sources. Currently the favoured 
explanation is in t erms of propagating acc retion rate fluctuations 
dLvubarskiil 1 19971 ; I Arevalo & Uttle^dl2006l) ■ although any similar 
scheme for multiplicative modul ation of random fl uctuations will 
produce a linear rms-flux relation dUttlev et al .120051) . In this model, 
random fluctuations in the viscosity occurring on all spatial scales 
in the accretion flow modulate the accretion rate further in, but high 
frequency variations occurring on large scales are damped as they 
move inwards. Irrespective of this, the presence of the linear rms- 
flux relation in more than one ULX suggests the same underlying 
variability mechanism is at work in ULXs as in the sub-Eddington 
XRB and AGN. 


5.2 Interband X-ray time delays 

The so ft lag in NGC 5408 X-1 was first detected bv iHeil & Vaughatd 
d2010l) . based mainly on the 2006 XMM-Newton observation of 
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NGC5408 X-1. IPe Marco et alj l l2013bh subsequently used six 
XMM-Newton observations (2006, 2008,2010a,b and2011a,b) and 
recovered similar soft X-ray lags in each, above frequencies of 
~few mHz. We extended the lag-frequency analysis to lower fre¬ 
quencies (~ 0.1 mHz) and find that the soft lag extends down to 
~ 1 mHz with an amplitude of tens of seconds, but neither the 
phase lag nor time lag is constant with frequency, with the phase 

lag constant at-0.2 rad (soft lag) above ~1 mHz but increasing 

at lower frequencies. A similar lag-frequency spectrum is found for 
NGC 6946 X-1. 

The soft vs. hard band coherence is high, indicating the vari¬ 
ations in each band are linearly correlated, and probably driven by 
the same process. Any additional source of uncorrelated variations 
adding to one band must be weak. The lags may then be imposed 
by the action of linear filters (“impulse response functions”) acting 
on the driving variations in one or both bands. 

The interpretation of these lags remains challenging. In the 
following, we discuss possible implications and comparison with 
other accreting sources. 


5.2.1 Comparison with AGN scaling relations 


Soft X-ray time lags are now well established in AGN (see Uttlev 
et al. l20l4 for a review), where they are observed tvpicall^at fre¬ 
quencies above ~ 0.1 mHz, below which a hard lag dominates. The 
most accepted scenario for AGN assumes that the hard lag is pro¬ 
duced by accretion rate fluctuations moving inwards through the 
accretio n disc and energising a radially extended corona (e.g.. Ko- 
tov et al. l200lU Arevalo & Uttlevll2006h . The soft lags are thought 
to be produced by a separate process - “reverberation”, i.e. light- 
travel time delays as the primary coronal emissi on is reprocessed in 
the inner regions of a partially ionise d disc (see lFabian et al.ll2009L 
and the review bv lUttlev et al.ll20l4 and references therein). Typi¬ 
cally the reverberation is thought to be generated by the inner disc, 
within ~ lOr,, of the c e ntral bl ack hole. 

De Marco et al. ( 12013ah studied a sample of 32 AGN and 
found a scaling relation between the black hole mass and the soft 
time lag using the 15 objects where a soft time delay was detected. 
If their scaling relation extends to ULXs, we would expect the soft 
lag to be located at frequencies in the range 20-200 mHz with am¬ 
plitudes of 0.03 - 0.4 s for black hole masses of 10 - IOOOMq, as 
expecte d for stellar-mass to inte rmediate-mass black holes (see e.g. 
Fig. 9 of iDe Marco et al.l2013^ . As shown in Fig.|3 the soft lag in 
NGC 5408 X-1 occurs at frequencies of the order of mHz with am¬ 
plitudes of tens of seconds, much larger than expected for a stellar 
mass black hole binary and more typical of a low-mass AGN. 

However, NGC 5408 X-1 is clearly identified with a stellar 
counterpart and associated nebulosity on the outskirts of the dwarf 
galax y NGC 5408 and is therefore highly unlikel y to be an AGN 
(e.g.. [Pakull & Mirionil [ 200 ^ ICseh et alj|2013h . Indeed, its X- 
ray spectrum and coarse timing properties fit very nearly into the 
range of kn own behaviours for ULXs in the proposed ‘ul tralumi- 
nous state’ jGladstone_et_alJ|2009yMiddletra^£t_ayj20]T|^Sutton 


et al. l2013l) , indicating it is likely to be a relatively sm all black hole 
accreting at super-Eddington rates dMotch et alj|2014l) . 


5.2.2 Association with XRB QPOs 

Another suggestion discussed bv IPe Marco et alJ ( l2013bh is that 
the soft lag is associated with a QPO, as sometimes occurs in BH 
XRBs. There are some problems with this interpretation. Foremost 


is that the QPOs claimed for NGC 5408 X-1 are generally not clear 
and distinct features like QPOs in BH XRBs (of the “type-C” va¬ 
riety). We will discuss this further in a companion paper (Vaughan 
et al., in prenJjbuteven_j|_^e_^akethe_^SD_firs_of^Dheeraj & 
Strohma^er(|2012l|jD^Marco_et_a^jj20]3b|j_o^Caballero-Garcla 
et al. ( I2OI3I) at face value, the soft lags extend over a much broader 
range of frequencies than the QPO-like features, indicating the lags 
are associated with the noise PSD. It is also not clear that lags in 
ULXs, e.g., between the 0.5-0.9 keV and 1.2-7 keV energy bands 
at ~1 mHz QPOs, can be compared to those in XRBs, which are 
usually observed between hard er energy bands (e.g., 3-5 keV and 
5-13 keV) for QPOs at ~1 Hz dRemillard et al.ll2002l:l^ella et al.l 

|20o4 . 


5.2.3 Possible origins for the soft lag 

The most popular explanati on of soft X-ray time lags in AGN is i n 
terms of reverberation (e.g.. lFabian et alj|2009l : lUttlev et alj|2014h . 
Can a similarexglanatmn_^oldjnJhese_^L2C^Caballero-Garcla 
& Fabian TioiOl) fit the XMM-Newton spectrum of NGC 5408 X-1 
using reflection models. Their best-fitting model includes reflec¬ 
tion from an ionised disc extending in to 2.5r^. Their model 
accounts for the “soft excess” in terms of emission from the disc 
(e.g., ionised N, O, Fe-L emission lines and recombination con¬ 
tinual, rather than using an additional continuum component, and 
so represents an upper limit on the strength of a reflection compo¬ 
nent. 

This model does not easily explain the timing properties. If the 
soft lag is the result of strong soft X-ray reflection responding to the 
primary X-ray luminosity, then the soft excess should be present in 
the covariance spectrum, but the observed covariance spectra for 
NGC 5408 X-1 and NGC 6946 X-1 both show no soft excess ( Mid- 
dleton et al. l2015h . Furthermore, the magnitude of the delay would 
imply an unreasonably large Mgu (see Sect. l5.2.Tt . Most of the re¬ 
verberation signal in the reflection model is produced in the inner 
< lOrj. If the observed maximum soft time lag of 10 s is not 
“diluted” (due to the contribution of direct emission to both bands), 
then the delay corresponds to a path length of ~ 3 x 10® km. This 
implies a distance of ~ 2 x 10''(MBa/100Mo)rj from the illuminat¬ 
ing X-ray source. Unless the black hole mass is > 10® Mg, com¬ 
parable to a low mass AGN, the reverberation must be originating 
at » lO^r-j, inconsistoi^_witlMhe_^eflectiogjgodel_^^Caballero- 
Garcla & Fabian T201(]|) . Since the effect of dilution is to make the 
reflection site even more distant in order to produce the observed 
delay, it would most likely produce strong emissions lines (e.g., Fe 
Ktr) that are relatively narrow, which are not seen in ULX spectra 
(e.g.. lCTadstone et al.l2009l : ISutton et alj20l3) . 

_ I ndeed, recent NuSTAR ob servations (e.g., iBachetti et al.l 

I 2 OI 3 I; IWalton et alJl20l412015bll4 have demonstrated that ULX 
spectra lack the obvious signatures of strong reflection (such as 
iron emission lines and the “Compton hump”) commonly seen in 
Sevf ert Is that display a soft lag reverberation signal (e.g.. Zoghb i 

et al. I2OI lUKara et al.ll2013ll . In particular. IWalton et alJ l l2015bl) 

studied NuSTAR, XMM-Newton and Suzaku observations of Holm- 
berg IIX-1, this ULX being the closest analogue to NGC 5408 X-1 
and NGC6946 jyjto^as^een_observedb 2 ^feS^^JCaballero- 
Garcia & Fabian l2010l : ISutton et alj|2013l) . and found no evidence 
of reflection in the spectrum. Hence a similar reverberation origin 
for the soft lags in ULXs appears unlikely, as reflection is not a 
strong component of the typical ULX spectra. We note that even if 
NuSTAR data would be useful to rule out (or not) the reverberation 
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origin of the soft lag, NGC5408 X-1 is too soft to make a good 
NuSTAR target. 

The spectrum of these “soft ultraluminous” sources such as 
NGC 5408 X-1 and NGC 6946 X-1 are often described in terms of 
two components_|Gladstone_e^_alJj2009ySutton_et_^]Jl201^Mid- 
dleton et al. 120151) . The softer component, dominating below ~ 1 
keV has a quasi-thermal spectrum and shows little, if any, short 
timescale variability. The latter fact is inferred from the lack of the 
“soft excess” in the rm s and covariance spectra (see Fig.[T]and also 
[Middleton et alj[2015h . The harder component, do minating over 
~ 1 - 10 keV , resembles a cut-off power law (see also lWalton et al.l 
I2OI3II2OI4I) , and produces most or all of the rapid variability. In 
the context of the supercritical accretion model (Shakura & Sun- 


yaev 

1973HKin2 et alJl200ll:lPoutanen et al J 20071 iMiddleton et alJ 

I2015 

) the hard component represents emission from the inner ac- 


cretion flow, while the soft component is thermal emission from the 
base of an optically thick, massive wind driven off the disc at larger 
radii. 


The lack of variability of the softer component suggests it 
plays no role in generating the soft lags. Further, the high coher¬ 
ence between soft and hard bands is most simply explained if there 
is a single “driver” of the variability. Together, these are consistent 
with the soft lag being intrinsic to the harder component and not 
a delay between the two spectral components. In other words, it is 
only the hard component that varies rapidly, and its variations on 
timescales shorter than ~ 1 ks occur at higher energies first, with 
the softer emission (from the same spectral component) taking up 
to ~ few seconds to respond. These delays may be intrinsic to the 
emission mechanism (from the inner accretion flow) or imposed by 
processes intercepting and delaying some fraction of this emission. 

One explanation for X-ray time lags is in terms of scattering in 
an intervening medium. A soft lag could be produced as hard X-ray 
photons from the primary X-ray source pass through and are down- 
scattered to lower energies. Low energy photons typically have un¬ 
dergone more scatterings and so escape the scattering medium af¬ 
ter a longer delay. The lag-frequency spectrum for NGC 5408 X-1 
shows clear frequency dependence, being constant in neither phase 
lag nor time lag. Such lags are however difficult to reproduce in 
a simple scattering scenario which would more naturally produce 
an approximately c onstant time lag (at frequencies below th e wrap 
aroun^_^g^enc^ Miller_et_jL|j^^^JZogtoi_e^jyj201^Uttley 


et al. l20l4 . 

A further possibility is that the soft lag somehow results from 
the propagation of photons through an optically thin shroud of ma¬ 
terial, likely to result from an expanded wind (a prediction of the 
super-critical ULX model). Should the absorption opacity of this 
material be low in the hard band and high in the soft band (per¬ 
haps as a result of high abun dances of Oxygen and Neon but rela¬ 
tively low abundances of Iron. fMiddleton et al.l20l4 . the hard pho¬ 
tons will arrive to the observer scattered (with the exact scattered 
fraction dependent on the Thompson optical depth) whilst the soft 
photons will be absorbed and re-emitted. There are indeed strong 
indications that residuals at soft energies in the time-averaged spec¬ 
tra are associated with absorption and emissio n features associated 
with a strong outflow l lMiddleton et al.l[20141 : Middleton et al. in 
prep) and, should this provide an origin for the soft lag, the magni¬ 
tude of the lag should provide constraints on the extent of the wind. 
However, such models incorporating outflows are by their very na¬ 
ture “messy” with the impact of reverberation and absorption heav¬ 
ily dependent on inclination angle and structure of the wind. 

A final possibility is that the lags are intrinsic to the X-ray 
source itself, i.e., that the part of the inner accretion flow that first 


responds to inward moving accretion rate variations has a harder 
spectrum than the later responding parts of the flow. This is in the 
oppo site sense to the models used to explain hard lags in XRBs 
(e.g.. lKotov et alj 20^ I Arevalo & Uttlevl2006l) . But we should per¬ 
haps not rule out this idea, as the structure of the hard X-ray emit¬ 
ting inner regions of ULXs may be quite different from the corona 
or jet-base thought to produce the thermal/non-thermal, variable 
hard X-ray spectrum in sub-Eddington XRBs and AGN. 


5.2.4 Low frequency lags 


At frequencies below ~ 1 mHz the sign of the time lag estimate 
changes and becomes a hard (or positive) lag. We do caution that, 
even with multiple observations of one of the brightest ULXs, the 
lag at ~ 0.1 mHz has large uncertainties, with the lag exceeding 
zero only at the 2 - 3cr level. 

Hard lags at low frequencies, along with soft lags at higher 
frequencies, have been observed in many AGN 1 Fabian_et_alj20^ 


Zoghbi et al.l201 lUErnmanoulopoulos et alj201 it Kara et alj2013l : 


De Marco et al.ll2013al) . Th e best studied source is tbe narrow lin e 

Seyfert 1 1H0707 -495 (e.g.. lFabian et al.l2009l:l^ghbi et aliTOUl) . 
[Kara et ^ ( 1201 3l) argued that the different energy dependence of 
the high and low frequency lags (soft and hard, respectively) in this 
object revealed that different emission processes are involved. In¬ 
terestingly, the lag-energy spectra for NGC 5408 X-1 at low and 
high frequencies (Eig. II look like very similar (after a change 
of sign and scaling). But given the low significance of the low 
frequency hard lag it would perhaps be premature to place much 
weight on this. 


The limited spectral-timing data for ULXs leaves few clues 
about the origin of the lags. Longer observations with higher S/N 
would be able to access both lower and higher frequencies, and 
allow us to better estimate the frequency and energy dependence of 
the lags. Such constraints will be crucial for ruling out lag models, 
but may require future, larger X-ray missions. 
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NGC 5408 X-1 (2010a,b) NGC 5408 X-1 (20na,b) 




NGC 5408 X-1 (2014) 



Figure Al. rms-fiux relation for NGC 5408 X-1 in the soft (0.5-0.9 keV) and hard (1.2-7 keV) energy bands, for observations from 2010, 2011, and 2014. 
The rms is measured over the 4-50 mHz frequency range from segments of length 150 s. The dashed lines show the 95% “confidence bands” around the 
best-fitting linear model. 


APPENDIX A: RMS-FLUX RELATION 

In this appendix we provide the plots of the rms-flux relation for NGC 5408 X-1 and NGC 6946 X-1 for each year separately. The results are 
presented in Sect.[3l 


APPENDIX B: CROSS-SPECTRUM 

In this appendix we present the cross-spectral analysis for NGC 6946 X-1, including the noise-subtracted coherence, phase lag, and time lag 
between the soft and hard energy bands. The trend of the lag is similar to that observed in NGC 5408 X-1, although the coherence is badly 
constrained, and both the soft and hard lag are below the 3cr level. Therefore, we cannot confirm the lags in this ULX. 
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NGC 6946 X-1 (2007a,b) 



Figure A2. rms-flux relation for NGC 6946 X-1 in the soft (0.5-0.9 keV) and hard (1.2-7 keV) energy bands, for observations from 2007. The rms is measured 
over the 4-50 mHz frequency range from segments of length 250 s. The dashed lines show the 95% “confidence bands” around the best-fitting linear model. 
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Figure Bl. Results of the cross-spectral analysis between the soft (0.5-0.9 keV) and hard (1.2-7.0 keV) energy bands for NGC 6946 X-1 showing, from upper 
to lower panels: the noise-subtracted coherence; phase-lag; and time lag (the inset panel shows a zoom-in of the soft lag). The middle panel shows a constant 
model (dot-dashed line) and a power law plus a constant model (dashed line) we fitted to the phase lag (see text). 





















































